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Abstract —We consider a single-cell massive MIMO full-duplex 
wireless communication system, where the base-station (BS) is 
equipped with a large number of antennas. We consider the 
setup where the single-antenna mobile users operate in half- 
duplex, while each antenna at the BS is capable of full-duplex 
transmissions, i.e., it can transmit and receive simultaneously 
using the same frequency spectrum. The fundamental challenge 
in this system is intra-cell inter-node interference, generated 
by the transmissions of uplink users to the receptions at the 
downlink users. The key operational challenge is estimating and 
aggregating inter-mobile channel estimates, which can potentially 
overwhelm any gains from full-duplex operation. 

In this work, we propose a scalable and distributed scheme to 
optimally manage the inter-node interference by utilizing a “one- 
hop information architecture”. In this architecture, the BS only 
needs to know the slgnal-to-interference-plus-noise ratio (SINK) 
from the downlink users. Each uplink user needs its own SINK, 
along with a weighted signal-plus-noise metric from its one-hop 
neighboring downlink users, which are the downlink users that it 
interferes with. The proposed one-hop information architecture 
does not require any network devices to comprehensively gather 
the vast inter-node interference channel knowledge, and hence 
significantly reduces the overhead. Based on the one-hop informa¬ 
tion architecture, we design a distributed power control algorithm 
and implement such architecture using overheard feedback 
information. We show that, in typical asymptotic regimes with 
many users and antennas, the proposed distributed power control 
scheme improves the overall network utility and reduces the 
transmission power of the uplink users. 

I. Introduction 

With each generation of wireless standards, the number of 
antennas at the infrastructure nodes has continued to grow 
to meet the increasing demand of mobile data. For example, 
both cellular and WiFi standards now support up to 8 antennas 
at infrastructure nodes. A promising and now standardized 
approach to use multiple antennas at the infrastructure is to 
use Multi-User Multiple Input Multiple Output (MU-MIMO) 
for supporting multiple uplink or downlink data streams in 
same time-slot. More recently, massive MIMO regime has 
been explored, which allows for a large number of antennas 
to reside at the infrastructure (with orders of magnitude 
more antennas compared to conventional MIMO systems), has 
attracted large interests in both academia and industry. The 
theoretical ( ifTSl . ifTSll . lfT9l ) and experimental results (123, 
E3) have shown that massive MIMO can lead to a number 
of desirable properties from the network design point of view, 
that include reduced inter-beam interference, improved energy 
efficiency, and reduced inter-cell interference, among others. 
As a result, massive MIMO is being considered for standard¬ 
ization as one of the key technologies in next-generation of 
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wireless systems in both cellular and WiFi (e.g., 3GPP LTE 
Release 12 |[2l and 802.11-ax |[3l). Channel models with 64 
BS antennas have already been standardized El, with research 
platform supporting even larger configurations; e.g. 96 BS 
antennas ll26ll . 

In addition to MU-MIMO, another potential avenue to 
achieve higher spectral efficiency is to leverage full-duplex 
transmission, where a full-duplex-capable device can transmit 
and receive at the same time using the same frequency 
spectrum. The full-duplex mode of network operation has 
the potential to double the spectrum efficiency of wireless 
networks and can bring substantial flexibility to higher layer 
design E). In fact, in-band full-duplex has already become part 
of the ongoing standard both in 3GPP HI and 802.11-ax |f3| . 
While it can be hard to integrate the full-duplex capabilities 
to client devices due to the processing and energy constraint, 
it is in fact feasible to design near-perfect full-duplex base- 
stations thanks to available freedom (bigger size, non-battery- 
powered operation) in their designs (e.g., see Ha, na, Ea 
and the references therein). One method to leverage full- 
duplex infrastructure with half-duplex mobile handsets is to 
have simultaneous multi-user up- and downlink transmissions. 
However, the potential for simultaneous up- and downlink 
MU-MIMO at the BS leads to a new challenge - the inter¬ 
node interference within each cell, i.e intra-cell interference, 
as the transmissions of uplink interfere with the receptions of 
downlink as illustrated in Figure HI The inter-node interference 
hence poses a fundamental challenge to enabling full-duplex 
at the BS and needs to be managed efficiently. 

In the presence of intra-cell interference, a centralized 
scheme can be used where the infrastructure device aggre¬ 
gates comprehensive knowledge of all the uplink, downlink 
and inter-node channel to make jointly optimal decisions on 
resource allocation. The centralized scheme may be useful for 
sparse scenario with small number of users. It will, however, 
incur a significant amount of overhead as the number of users 
grows, with the inter-node interference being the dominant 
bottleneck. For example consider the setup with M full- 
duplex-capable antenna^j at the BS, and with Kui half-duplex 
single antenna uplink and Kdi half-duplex single antenna 
downlink users, henceforth denoted as (M, Kui, Kdi) massive- 
MIMO full-duplex system, in FigureHl The centralized control 
scheme requires the full-duplex BS to aggregate knowledge 
of Kui X Kdi inter-node interference channels, in addition to 
M X Kui uplink and M x Kdi downlink channels. Note that 
there has been recent development on uplink/downlink channel 
estimation in Massive MIMO with reasonable complexity HTll . 

’a full duplex capable antenna is capable of transmitting and receiving at 
the same time. It can be realized by either using one antenna (e.g., H) or 
using a pair of antennas (e.g., (mnu) 


EH. For example, the M x Kui uplink channels can be ( 
mated by letting K^i uplink users send pilots to the BS, and 
downlink channels can be obtained by letting Kdi down 
users send pilots and estimated by utilizing channel recipro^ 
Hence the M x (Kui + Kdi) up- and downlink channels cai 
obtained with {Kui + Kdi) transmissions of overhead (pile 
However, to obtain K^i x Kdi inter-node interference chant 
each of the Kui x Kdi channel information needs to estim; 
at the uplink/downlink users and sent to the infrastructure; 
inter-node interference channels are two-hop away from 
BS. Hence acquiring all the channel knowledge of Kui x 
inter-node interference channels incurs overhead at leas 
the order of Kui x Kdi- As an example, under the seen 
with Kui = 15 uplink and Kdi = 20 downlink users, tl 
can be as many as Kui x Kdi = 300 inter-node interference 
channels, which is significantly more than Kui + Kdi = 35. 
Therefore, for centralized scheme, obtaining the inter-node 
interference becomes the main bottleneck and will render the 
overall centralized architecture non-scalable. 

In this paper, we thus seek distributed methods for resource 
allocation that significantly reduce who needs to know what 
about which channels. We propose, to the best of our knowl¬ 
edge, the first scalable and distributed interference manage¬ 
ment techniques in massive MIMO full-duplex systems with 
half-duplex clients. Our contributions are listed as follows. 

• We create a scalable information architecture, namely 
the one-hop information architecture, that specifies the 
network information needed in each cellular devices to 
perform optimal interference management. Specifically, 
we consider an architecture where (i) BS knows only 
SINR of downlink channels, much like current cellu¬ 
lar systems which rely on periodically collected SINR 
information for power and rate control, and (ii) up¬ 
link transmitters know only about interference channel 
strength, along with a weighted interference-plus-noise 
(IN) metric, from the nearby downlink receivers. There¬ 
fore, for the {M, Kui, Kdi) massive-MIMO full-duplex 
system, compared to a fully centralized decision, the 
proposed system requires only local, at most Kdi one- 
hop information at each uplink user from its neighboring 
downlink receivers, thus avoiding the significant non- 
scalable overhead for aggregating the Kui x Kdi inter¬ 
ference channel information at the BS. 

• Based on the one-hop information architecture, we pro¬ 
pose a distributed power allocation algorithm to achieve 
the optimal sum-utility in the system. By analyzing 
the structural properties of the {M, Kui, Kdi) massive- 
MIMO full-duplex system, we identify that the up- and 
downlink power control, which is coupled by inter-node 
interference, can in fact be carefully decoupled. Such a 
decoupled structure turns out to be an enabler for the one- 
hop information architecture and the optimal distributed 
algorithm. 

• We further show that, in typical asymptotic regimes where 
the amount of antennas scales up with the number of 
users, distributed power allocation can not only improve 
the overall system utility, it can also asymptotically 
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Fig. 1. Full-duplex multi-user network with Kui uplink users transmitting 
to the BS and K^i downlink users receiving from BS. 


reduce the overall amount of transmission power and 
hence improve the energy efficiency of mobile users. 

Note that, there can be scenarios with dense users and hence 
strong inter-node interference, where half duplex mode can 
be more efficient solution. The scope of this paper is not to 
specify when and where to use full duplex mode of operation, 
but rather to explore the optimal and scalable power control 
schemes for massive MIMO when the full duplex mode is in 
use. 

The rest of the paper is organized as follows. In Section Hill 
we formulate the physical layer model of massive MIMO 
full-duplex system. We proceed to present the power control 
problem in Section |IVl In Section lYl we present the one- 
hop information architecture that facilitates us to propose dis¬ 
tributed interference management algorithm. An overhearing- 
based scheme is presented in Section [Vl] to implement the 
information architecture. We present asymptotic analysis in 
large number of antenna regime in Section IVIII Numerical 
evaluation are provided in Section rVIIIl followed by conclusion 
remarks in Section HXl 

H. Related Works 

There has been extensive research on resource allocation in 
various interference-constrained wireless networks. The inter¬ 
ference in the wireless networks is often captured by the SINR, 
which in turn determines the achievable transmission rate. 
Power control is an important mechanism to mitigate network 
interference and have been well studied (see survey CD- 
Such problem often requires solving nonconvex optimization 
problem that is proven to be NP hard to solve CD. 

One line of research has focused on centralized control 
where a centralized controller gathers all the information 
across the network to perform resource allocation and in¬ 
terference management. For instance, the authors in ca 
proposed centralized power control algorithm under high SINR 
regime using geographic programming. Signomial program¬ 
ming based approaches are proposed in El to achieve a local 
optimal solutions. An iterative optimal algorithm was proposed 
in 1201 to achieve global optimality for a subset of non¬ 
convex objective functions. A branch and bound algorithm 
was proposed in l29l to solve weighted sum-rate maximization 
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problems. Approximation algorithm was proposed in Il28l to 
solve the NP-hard power control problem. 

Another line of research has studied distributed power 
control algorithms for interference management. In the high 
SINR regime, distributed power control algorithm was con¬ 
verted to geometric programming 0. Distributed price-based 
power control algorithm was proposed in Cl and shown to 
be optimal in certain scenarios. Distributed Power Control 
and scheduling based on back-pressure algorithm is proposed 
in in high SINR regime. Recently, optimal distributed 
scheduling based on SINR model is proposed in to achieve 
transmission rate over a power control region 0. 

Some recent work studied resource allocation for full du¬ 
plex systems. The authors in ifTol studied optimal subcarrier 
allocation for Full Duplex OFDMA network. A cell parti¬ 
tioning method was proposed in ll24ll that allocates frequency 
resources according to the partitioning of a single cell. 

III. System Model 

In this section, we describe the channel model for a single¬ 
cell massive MIMO full-duplex system where the BS is full- 
duplex capable while the mobile users operates in single 
antenna half-duplex mode. We assume the full-duplex BS is 
equipped with M antennas to serve both uplink and downlink 
users in the same time/frequency slot. The sets of uplink and 
downlink users are denoted by /C„; and K-di respectively, with 
|/C«/| = Kui and \lCdi\ = Kdu as shown in Figure □ 

For each uplink/downlink user k, an associated utility 
function Uk{r) represents the utility it achieves under trans¬ 
mission rate r. It reflects the level of satisfaction achieved 
for each of the users at a given rate. Each utility function is 
assumed to be a non-decreasing, continuously differentiable, 
strictly concave function, where is the derivative of the 

utility function. The concavity property captures the dimin¬ 
ishing return property of the utility achieved by the users 
in the network ll^ . These utility functions are also used 
to capture different fairness requirement over the network. 
We are interested in controlling the transmission power, over 
both uplink and downlink in the massive MIMO full-duplex 
system, so that the achievable transmission rate provides the 
maximum overall sum-utility over the network. To characterize 
the physical-layer relationship between the transmission power 
and the corresponding achievable transmission rate, we start 
by presenting the channel model of the massive MIMO full- 
duplex system. 

The propagation channel model in our system consists of 
two parts; small-scale fading due to multipath and large-scale 
fading due to path loss. The channels are represented by the 
set Q = {G“*, G^}, where each element in Q denotes 

the propagation matrix of complex-valued channel coefficients 
for uplink, downlink and inter-node interference channels, 
respectively. We can express G^, k € as 

Gfe = Lfc(Dfc)^ 

where represents the small-scale fading matrix with a 
magnitude of one for each entry, and (Dg )2 = diag{\f^) is 
a diagonal matrix whose diagonal entries constitutes a vector 
g^' of path loss coefficients. 


A. Uplink 

At the BS with M antennas, the received uplink signal y“* 
at the BS is an M X 1 vector given as 

= G“'x“'(1) 

where G“^ is an M x K^i uplink channel matrix and 
is a Kui X 1 vector, with each element representing the 
transmit signal of the uplink users. The vector z“^ represents 
the receiver additive noise, each entry is independent and 
identically distributed (i.i.d.) and drawn from a circularly 
symmetric complex Gaussian distribution with zero mean and 
variance of Nq denoted as CAf{0, Nq). 

In the uplink, we let P“^ = , • • • , ] denote 

the vector of transmit powers, where the represents the 
transmit power of uplink user i. Each user is subject to a 
maximum instantaneous power constraint i.e., < 

,V* e /C„,. 

Eollowing the techniques in IfTSlI , 1211 . the BS can simply 
process the received signal using maximum-ratio combin¬ 
ing (MRC), 

(G“')^y“' = ^ (2) 


where superscript H denotes conjugate transpose. We hence 
have 



where (D“^)^ is a K^i x K^i matrix and is an M x Kui 
matrix. In massive MIMO, the number of antennas at BS is 
much larger than that at mobiles. Thus when M ^ K^i, under 
favorable propagation conditions as manifested in ITSl . the 
columns of the small-scale fading matrix are asymptotically 
orthogonal, the following convergence is obtained ITSl . 
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Now substituting (HI into (|2l), we have 
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We then compute the Signal-to-Interferene-and-Noise ratio 
(SINR) for uplink user i as 


SINRf 


MPfgf 

No 


( 6 ) 


where a/ gf’- is the ith entry in the diagonal matrix (Dg*)5. 
Note that in the above expression, we have used the asymptotic 
result in (|5]l. Under the domain of massive MIMO, we assume 
the SINR expression is valid for large value of M, hence 
ignoring the negligible error term. 

Therefore, with receiver MRC, the uplink user i achieves a 
transmission rate 


'Ul ^ul s 
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^ul ; 


(7) 


where, recall that, is the corresponding path loss coefficient 
for i-th user. 









4 


B. Downlink 

The received signal at the Kdi downlink users is a 
combination of the downlink signals and the interfering uplink 
signals, and is given by 


denote 


where y“' is a Kdi x 1 received signal vector, is a 
Kdi X M downlink channel matrix and is an M x 1 
downlink transmit signal vector; is a Kdi x K^i inter¬ 
node interference channel matrix; is the receiver additive 

Gaussian noise which contains i.i.d. CM{0,Nq) entries. 

In the downlink, we let , • 

the vector where the element P^^ represents the transmit 
power to downlink user j. The downlink transmission power 
allocation is subject to a constraint on the total amount of 
transmission powers, i.e., = Ptot- 

Note that we assume perfect serf-interference cancellation 
at the full-duplex-capable BS. This is inline with recent results 
that self-interference cancellation is capable of suppressing 
self-interference below the noise floor da. 

Similar to the uplink, in massive MIMO downlink, under 
the conditions of M ^ Kdi, we have 


Gdlf^Qdl^H 

W 


D 


dl 

g ■ 


(9) 


M-3>Kdi 

Following the techniques in f2T\ . BS can perform conjugate 
beamfoming precoding, 

1 






( 10 ) 


where Dp = diag{'/P^) is a Kdi x Kdi diagonal matrix. The 
Kdi X 1 vector contains downlink messages where each 
message has unit power, i.e., E{\sf\) = 1 for all j G ICdi- 
We substitute (|9]l and (flOl i into ® and we obtain that 


V ^/M ) m:»k^i 




,dl 


G^x“' 

VM y/M' 


( 11 ) 


Because the uplink and downlink transmissions are over the 
same frequency band, the uplink transmissions causes interfer¬ 
ence to the downlink receptions, as shown in Figure [T] The 
interference channel strength between uplink i and downlink 
j is denoted as gfj. For the uplink user i, we let Mi C 
denote the set of downlink users that receives non-negligible 
interference from user i. For downlink user j, we let Mj C 
denote the set of uplink users that imposes non-negligible 
interference to user j. 

In the massive MIMO scenario where M ^ Kdi, the SI NR 
of downlink user j can then be calculated as 


SINRf = 


M ■ Pfgf 


E nl P 


ml 


-Nn 


( 12 ) 


where ygf is the jth entry in the diagonal matrix (Dg^js, 
and, similar to the uplink, we have ignored the negligible error 
term in (O. 

We consider the setup where the downlink receivers treat 
the inter-node interference as noise. Therefore, with transmit 


conjugate beamforming, the downlink transmission rate to user 
j € ICdi is expressed as 


which depends on both uplink and downlink power allocation 
schemes p“^ We henceforth denote the Interference-plus- 
Noise (IN) by IN, = 


IV. Power Control Problem 

In this section, we formally define the power control 
problem. The objective of power allocation is to control the 
transmission powers p^^^p^^^ to jointly manage the inter¬ 
node interference and intelligently allocate the downlink trans¬ 
mission powers under the interference, so that the overall 
sum-utility in the system is maximized. Recall that 
and Jidi-(^pdi^pui'^ defined in 0 and ( fT3] l respectively. 

The power allocation problem is defined as the following 
optimization problem. 

^max^^ ^ M{Rf{Pr‘))+ Y. U,{Rf{Pf,V-^)) 

i^fCui j^Kdi 

(14) 

s.t., Y Pf ^ PtL (15) 

j£lCdi 

0 < , g I^^^^pdi > 0^^. g (16) 

As previously discussed, the approach of centralized con¬ 
trol can be used to solve the problem at the infrastructure, 
where the full-duplex BS aggregates comprehensive knowl¬ 
edge across the cellular network. To optimally manage the 
inter-node interference imposed by the uplink users to down¬ 
link receivers, the infrastructure needs the knowledge of all 
the inter-node interference channel gains. From the form of 
optimization problem, the required amount of information at 
the BS is listed in Table |T] 0 


Information at BS 

Amount 

Uplink path loss: 

All i G K,,i 

Downlink path loss: g^'' 

All j G K-di 

Inter-node interference: 

All (i,j) pairs, iG/Cui,jGA/i 


TABLE I 

Required information for centralized algorithm at BS 


As discussed in the introduction, the above centralized 
solution is non-scalable primarily because of the significant 
overhead to potentially collect K^i x Kdi inter-node inter¬ 
ference channels. Therefore, our goal is to design a scalable 
distributed structure, where the network-wide information are 
not aggregated at the BS. Instead, the information is distributed 
across the cellular network so that each device only acquires 
information at most one hop away. 

^Note that in the above table, we do not include the parameters M, 
PfatPrnax’ which are assumed to be predetermined information that does 
not need to be communicated network wide. 
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V. Distributed Solution using one-hop 

INEORMATION 

In this section, we present a scalable network state infor¬ 
mation architecture, denoted by one-hop information. Such 
information architecture significantly reduces the overhead 
associated with centrally aggregating all the information at the 
infrastructure. Specifically, the one-hop information architec¬ 
ture is described below in Table HJlInl where the Interference- 
plus-Noise (IN) metric, i.e., rriij [f] is formally defined in (l29l l 
in Section IV-CI 


Information at BS 

Amount 

Downlink SINR SINRf [t] 

All j e Kdi 


TABLE II 

One-hop information architecture, BS 


users in the neighborhood in the range of the cell. Corre¬ 
spondingly, the high SINR model gives the transmission rate 
expression 


M■Pfgf 


=log 


( SiGAf, dijPi"'' + No ^ 


, 3 S Kdu 


( 18 ) 


which, similar to the uplink, is valid for Pj** > Pj. The 
parameter Pj is a small value as the number of BS antennas 
scales up. 

Following the techniques in lb), [SOl, we perform a change 
of variables by letting P“^ = log(P“^) and P^^ = log(Pj**) 
for all i G K-ui and j G K-di- With the new set of variables 
and under high SINR approximation (fTTl i lfTsT l. the utility 
maximization problem (fT4l) - (fT6l) becomes 


Information at uplink user i 

Amount 

Uplink SINR: SINRf [t] 

I 

IN metric: rriij 1^] 

All j G Ni 


TABLE III 

One-hop information architecture, uplink user i 


It can be observed from the above tables that the BS requires 
only SINR from the downlink users, which does not incur 
additional overhead to current cellular systems where BS peri¬ 
odically collects SINR from downlink users for power and rate 
control. The BS hence does not need to know information two 
hops away, e.g., inter-node interference. For each uplink user 
i, it only requires the one-hop information, i.e., the IN metric 
rriij for all users j G Afi- It does not require information more 
than one-hop away, e.g., migj[t] for all users io ^i,j G . 
Therefore, the multiplicative overhead associated centralized 
solutions disappears. In the the following sections, we show 
that under massive MIMO systems with high power gain and 
hence high SINR, the above one-hop information architecture 
enables design of optimal distributed power allocation scheme 
to manage inter-node interference. 


max E E (19) 

iGlCui j^^di 

s.t., n - log(Mp“')-P“' < - logiVo, i G K-ui (20) 

r,-log(M5f)-Pf+log 4exp(P“')+iVo) < 0, 

ieMj 

j G ICdi (21) 

^ exp(Pf') < P,l' (22) 

jSlCdi 

log P° < Pf < log P“L P G (23) 

Pf > log P°, J G /Cd/, rfe > 0, fc G Kui U Kdi (24) 


where r = [rk]k(^K^iuKdi- 

Observatoin 1: The high SINR provided by massive MIMO 
system enables us to transform the non-convex problem (fT4l i- 
(fTbl) to the above convex optimization problem. 


B. Leveraging Full-Duplex Architecture to 
Reduce Overhead 


A. High SINR Model 

In massive MIMO full-duplex system, as illustrated in the 
rate expressions © and (fTST l. large power gain is brought to 
each user owing to the large number of antennas at the BS. The 
high power gain in massive MIMO brings high SINR values at 
each receivers. Under the high SINR, the uplink transmission 
rates hence can be approximated by 

~ / M ■ ul \ 

Pf (^“')=log(SINRf)=log , i G ICuG 

(17) 

for P“* > P° where the parameter Pq is the power allocation 
constraint such that the high SINR approximation is valid, i.e., 
SINRf » 1. Under massive MIMO setup, P^,i G are 
small values. As M scales further up, the value of P° scales 
down. 

In the downlink, the SINR expression given in (fT2li is high 
thanks to the large value of M compared to the number of 


One of the main characteristics and bottleneck of the mas¬ 
sive MIMO full-duplex system is the inter-node interference. 
Such interference not only can incur significant overhead in the 
centralized solution, but also couples up- and downlink power 
control decisions. Because the uplink transmission power 
causes interference to downlink reception, the downlink power 
allocation will also be affected accordingly. In this section, 
we present an interesting observation that there exists a way 
to decouple the interplay between up- and downlink power 
allocations, thanks to the uni-directional uplink to downlink 
interference in full-duplex system. This structure helps reduce 
the overhead associated with obtaining network information 
to achieve the optimal utility for the aforementioned prob¬ 
lem (Ug-dSl). 

For each uplink user i G K-ui, we associate a Lagrangian 
multiplier gf that corresponds to the constraint ( |20] |. A 
Lagrangian multiplier is associated with downlink user 
j G ICdi that corresponds to constraint (l2Tl i. The following 
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dual function is then obtained 


C. Power Control using one-hop Information 


^(q“ , q ) = max B(r, q“‘, q'*') + max V^i (P“', q“', q”^') 

r p„i 

+ max Vdi V qf log No. (25) 

pdl 


where 

i3(r,q“U"')=E {U.(r.)-9f-n)+j2 {Uj{r,)-qf-r,), 

(26) 

K/(P“',q“',q"')= ^ qf-logiMgf)+ ^ qfPf^ 

- X! ■ log ( X! dij exp(-^“') + A^o), 

j&ICdi *ga(, 

y,,(P‘'',q‘'')= ^ qf-\ogiM-gf)+ ^ qfpf, 

j^f^dl j^^dl 

with vectors q“' = [qf]ieK^i, = [qf]jeKdn P“' = 

Hence, given q“*, q"^*, the dual problem is decomposed into 
three sub-problems, which leads to the following immediate 
observation. 

Observation 2. The uplink power control sub-problem 
maxgui 14 ^(P“^ q“^ q'^O’ the downlink power control 
sub-problem maxg^, t 4 ^(P‘^^q ‘^0 completely decoupled 

as separate problems. 

Remark: The Lagrangian approach uses a similar technique 
compared with other works that considers power control prob¬ 
lems in general networks (e.g., JD). However, the full-duplex 
capability at the BS provides a unique decoupling structures 
for up- and downlink power control. In the full-duplex systems 
that we have considered, inter-node interference limits system 
performance, and up- and downlink power controls are coupled 
because of inter-node interference. The decoupling structure, 
however, interestingly reveals that these two power control 
schemes can actually be performed separately, given some 
pricing information. Such separation structure enables us to 
significantly reduce the information passing overhead to design 
the optimal power control policies. 

For the uplink and downlink sub-problems, we obtain the 
gradient expressions as, 

dPf = 

The sub-problem maxr i3(r, q“( q'^*) is the primal rate 
adaptation, which can be solved with r* = r* = 

where i G ICuuJ G ICdi- Hence the Lagrangian 
multipliers can be considered as ‘price’ for achieving a unit 
of up- and downlink transmission rate. 

We next design a distributed power control algorithm based 
on the aforementioned intuitions. 


We formally present the distributed power control and 
interference management algorithm next, where the notation 
[a :]5 stands for projecting the value x onto set S, and the IN 
metric is 


mij[f\ 


IN,W 


(29) 


Algorithm Distributed Algorithm for Massive MIMO Full- 
duplex Systems based on One-Hop Information 

1 : Initialization phase. At f = 0, set g“*[0] = 0, qf[0] = 0, 
Pf^[0] = \ogP-i,, Pf[0] = \og{P,t/Kdi), IN,[0] = No 
for all i G Kui, j G ICdi- Go to step 5. 

2 : Downlink power update at the BS. The BS updates the 
transmission power to the downlink users according to the 
gradient (l28l l. 


Pfft] 

Pf[t-l]+yqf[t-l] 


. l-Pf ,Pf lt-l]>P^ 


where 7 is the step size. 

3: Uplink power update at the users. The uplink users 
update the transmission power to the BS according to the 
gradient (IZTT i. 


Pf^[t]=[pf^[t-1]+^ . [qf[t-l] 
j&Ni 


log Po<Pp [t-l]<log PP' 

(30) 


4: Price update. The prices G ICdi are updated as 

follows. 



qflt - 1 ] + 7 [r, [f - 1 ] - log (SINR,[t - 

(31)' 


where rj[t] = (f/j) The uplink user i G ICui 

update its price by 


qf[t]= i]+7h[^-i]-iog(SiNRi[f-i])] 


g“‘[t]> 0 ’ 

(32) 


where r,[t] = {U') ^{qf[t]). 
5: t ^ t + 1. Go to Step 2. 


In the algorithm, the up- and downlink power is updated 
according to the gradient (l27l i (l28T l in item 2 and 3. The 
price update in item 4 is along the direction of gradient of 
dual function i9(q“*,q"^*). Also note that in uplink control 
algorithm (l30l l in item 3, the term 

my [f - 1 ] exp (^Pf‘[t-l]J = qf[t - 

denotes the contribution of inter-node interference from uplink 
user i to the overall interference experienced at downlink 
user j, weighted by the price of user j. Therefore if user i 
contributes significant interference to downlink user j whose 
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achievable rate is already very ‘pricey’, user i’s power will 
tend to decrease. 

Note that the above power control algorithm is implemented 
distributively at the uplink and downlink users. The distributed 
nature of the algorithm is facilitated by the decoupling prop¬ 
erty between up- and downlinks previously discussed, which is 
captured by the one-hop information. Specifically, because of 
the decoupling, given the price downlink power control 
is implemented independently from uplink users. Hence BS 
does not require uplink power information. Uplink user power 
control is implemented autonomously without awareness of 
the downlink power control policy. Note that the the one-hop 
information takes advantage of such structure so that all the 
network devices has sufficient information to perform optimal 
power control. 

The following theorem establishes the convergence of the 
algorithm. The proof follows the similar lines as 0 and is 
hence neglected. 

Theorem 1. Let (jf be the set of that optimally 

solves the problem (I23-(lll- For small enough 7 , we have 
limt^oo(P“'[f],P‘''[f]) e 

It also follows from 0 that the algorithm converges to 
optimum geometrically, given in the following theorem, where 
G and 

e[k]=\ Y, Y U.{Rf{Pf[k],P^'[k])) 

i^Kui j^lCdi 

- Y u,{Rf{pY))+ Y u,iRfipY,p^'n) • 

i&ICul j&tCdi 

Theorem 2. There exist constants c and j3 such that for all 
k, e[k] < c/3^. 


VI. Power Control Implementation 

USING OVERHEARING-BASED SCHEME 


We have previously seen how the one-hop information 
architecture facilitates design of distributed algorithm for 
managing inter-node interference. Next, we present how the 
architecture can be implemented using an overhearing-based 
scheme. 

At the end of each slot t, a feedback message is transmitted 
from each downlink user j G ICm to the BS. The feedback 
message fbj[f] from the j-th downlink user contains the 
following value. 


gfM gf[f].SINRf[f] 
\Nj[t] MPf[t]gf 


(33) 


where INj[f] is the interference-plus-noise measured at the j- 
th downlink user. Note also that downlink user j has all the 
knowledge to keep a copy of price qj'‘[t] according to (l32l i. 

The feedback message fbj[f] is received at the BS and is 
overheard by all the uplink users in the neighborhood Nj of j- 
th downlink user. We assume these feedback packets are sent 
with low rate and can be received collision-free at BS and 
users via orthogonal frequency/time resources. The BS and 
the uplink users then process the packets as follows. 


• From the embedded pilot symbols in the feedback mes¬ 
sage fbj[f], the BS estimates the downlink path loss The 
BS then uses fbj[t] to calculate SINR^^[f] at the downlink 
receiver j at time t, i.e.. 


SINRf [f] 


fb j[t]MPf[t]gf 
qf[t] 


(34) 


• The uplink users in the neighborhood Afj of downlink user 
j overhear the feedback message. From the embedded pilot 
symbols in the feedback message fbj [f], the uplink transmitter 
i G ffj estimates the interference channel gain gk associated 
with downlink receiver j G ICdi ■ The uplink transmitters i G 
ffj then, from ( l29t and (l33T l, use the decoded feedback bbj [f] 
to calculate the metric mijf] = fb^ [fJgT and update the 
power allocation according to ( |30] |. 

Therefore, with the overhearing-based feedback scheme, the 
BS obtains the required SINK information in Table 1, and 
each uplink user obtains the IN metric in Table 2. Note that 
in Table 2, uplink SINR can be obtained from 0 where the 
uplink power is known to the user and path loss coefficient 
can be estimated via BS broadcasting pilots at initialization to 
all users. 


VII. Asymptotic Performance of 
Power Control 

Massive MIMO full-duplex system brings a large amount of 
possible degree of freedoms, hence facilitates the capability to 
simultaneously serve a large number of uplink and downlink 
users. In this section, we study the performance gain associated 
with power control, as well as the asymptotic scaling of the 
optimal power allocation in massive MIMO full-duplex system 
when both the number of antennas and the number of users 
scale up. 

We hence consider a sequence of scenarios with expanding 
number of users and antennas, indexed by I with I G Z+. 
For the k-th scenario, the set of uplink users, downlink users, 
and the number of BS antennas are denoted as /C^j, /C^;, Mfc, 
respectively. 

To study the system behavior when more and more users 
are added to the system, we study the monotonic scenario 
where Mi < Mi+i, and the uplink and downlink set of users 
monotonically expand across scenarios, i.e., C C 
ICh C • • • and IC^i C IC^ C K.^ C ■■■. We hence have 
|/C^;| —>■ 00 , |/C^;| —00 and M; —>• 00 as / —>• 00 . 

We assumes the inter-node interference g(j,i G K-nuj S 
K-di follow certain distribution with expectation Wj[gP] = £ 
uniformly across (i, j) pairs. This assumption is valid for the 
typical scenarios where users randomly appear in geographical 
locations of a cell, or where inter-node interference arrives via 
random paths. 

We assume the following conditions hold true for the utility 
functions of the users in the network: 

(Al) For each user k, limr_>.oo U'i.{r) = 0, k £ ICui U ICdi- 

(A2) The utility function Ui of i-th user satisfies Ui G U, 
where U = {[/^(r), • • • , [/^(r)} is a finite collection 

of utility functions. 






Assumption (Al) is motivated by the diminishing return 
idea of the utility functions, and holds true for the well-known 
a-fairness, e.g, 

j,l-a 

l/i(r) = UJ^- -,a > 0,a ^ l,uj^ > 0, 

1 — a 

which includes many well-known fairness criteria such as 
proportional fairness, minimum potential delay fairness and 
max-min fairness as special cases that corresponds to different 
values of a. Assumption (A2) states that there are finite 
choices of utility functions for each user. These choices 
correspond to, for instance, different types of data traffic, 
quality of service requirement, etc. 

For ease of exposition, we henceforth use the following 
notations: 

• Vector P“*>*’* and respectively denote the optimal 

uplink and downlink power allocation that maximizes the 
network sum-utility (i.e., (fT9]l- (l24l i) in the l-th scenario. 

• The utility i € }C[^i U respectively represents 
the up- and downlink utility under optimal power allocation 

puz.i.* jjj jjjg scenario. 

• The utility and D: rep¬ 

resent the utility when each uplink user transmits at maximum 
power P^ax’ ™rl when the downlink users transmits at optimal 
power that maximizes the sum downlink utility under 

The next theorem states that when the number Mi of the BS 
antennas scales at the same order of the multiplication of |/C^; | 
and |/C^;|, and when the number of downlink users scales at 
the same order or faster than the number of uplink users, then 
asymptotically the optimal uplink transmission power scales 
down. 

Theorem 3. Suppose limj^oo . ; , = C for some 

II 

constant C >> 1, and limsup,. < oo. For any 

I'Vdll 

0 < p < 1, we let Qp C be the subset of uplink users 
with < pP^ax- Then limz^oo = 1- 

Proof: See Appendix A. ■ 

Remark: 

(1) Theorem [3 studies the asymptotic regime when 

lim;_>oo I , = C, where the BS antennas Mi scales 

at the same order of the multiplication of and AT^;. 
Intuitively this means the amount of antennas grows suf¬ 
ficiently fast to combat simultaneously the growing inter¬ 
node interference (from the growing number of uplink users) 
and also to provide sufficient degrees of freedom to the 
downlink users to guarantee a non-trivial SINR gain at the 
downlink users expressed in (fT2l l. Theorem [3] also assumes 
that limsup, -Lti < oo, i.e., in the asymptotic regime, the 
number of downlink users scales at the same order or faster 
than the number of uplink users is valid in typical scenarios 
where the network is often dominated by the downlink traffic. 

(2) Theorem [3 states that in the considered asymptotic 
regime, the optimal transmission power asymptotically scales 
down. This is intuitive because the uplink SINR grows lin¬ 
early with the number of antennas. However, because of the 


(a) 


(b) 



Inter-node Interference Channel Strength 
(C) 




Inter-node Interference Channel Strength 

(d) 


Inter-node Interference Channel Strength 



Fig. 2. Perfomiance evaluation of power control algorithm for two-user case, 
(a) Optimal power scaling with weighted proportional fair utility; (b) Optimal 
power scaling with weighted minimum potential delay utility;(c) Comparison 
with Naive scheme with = Pffax under weighted proportional fair 

utility; (d) Comparison with Naive scheme with under weighted 

minimum potential delay utility. 


diminishing decay property of the utility function, the marginal 
utility gain due to the increasing number of antennas will 
reduce. Therefore, decreasing the uplink transmission power 
(e.g., decrease power by p) will not bring significant loss in 
the uplink utility, it does, however, bring significant downlink 
utility gain. This is because, as the number of uplink users 
increases, the aggregated interference can be large at each of 
the downlink users. 

VIII. Numerical Evaluations 

We examine the power control algorithm performance for 
the massive MIMO full-duplex with one uplink user and one 
downlink user. We consider the scenario with 128 antennas at 
the BS, i.e., M = 128. The uplink and downlink path loss pa¬ 
rameters are = —60 dB, = —70 dB respectively. The 
power control is subject to the maximum power constraints 
Pmax = 23 dBm, = 45 dBm. The noise power is 
Nq = —30 dBW. These parameters are set in accordance 
with practical values in existing LTE standards. Eig. |2a) and 
Eig. |2c) correspond to weighted proportional fair utility func¬ 
tion, with = log(r), U^\r) = 2 log(r). Eigure|2b) and 

Eigure I2d) correspond to weighted minimum potential delay 
utility function, i.e., (7“*(r) = —l/r, = —2/r. Eor this 

scenario, it is optimal for the downlink user to always transmit 
at the maximum power since there is only one downlink user. 

Eigure|3!a)-(b) plot the variation of the optimal uplink trans¬ 
mission power with the gain of the inter-node interference 
channel. Eigure|3!c)-(d) compare the sum-utility under optimal 
power control and when = Pmax- observed in the 
figures, when the inter-node interference channel is weak, it 
is optimal to transmit at the maximum uplink power since 
the impact on downlink transmission is insignificant. Hence 
in this regime the uplink transmits power at the maximum 
power and there is a negligible loss at the downlink as evident 
in Eigure I2c)-(d). As the inter-node interference channel 
gain grows stronger, for both cases of utility functions, the 
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Fig. 3. Convergence of power control algorithm, (a) Proportional fair utility; 
(b) Minimum potential delay utility. 


optimal uplink transmission power decays. This is because 
under strong interference channel, transmission of uplink users 
imposes significant interference to the downlink transmission. 
Because of the diminishing return property of utility functions, 
reducing the uplink transmission power can alleviate downlink 
utility loss and bring higher overall sum utility in the network. 

Figure |2jc)-(d) also highlight the importance of power 
control in moderate to high interference regime. As the gain in 
the inter-node interference channel becomes stronger, power 
control is increasingly more important to manage the interfer¬ 
ence created by the uplink users to the downlink receivers to 
maintain the overall system-level utility, as illustrated by the 
increasing utility gains in Figure |2c)(d). 

Fig |2ta)-(b) plot convergence performance of the pro¬ 
posed algorithm. We consider the scenario with 128 an¬ 
tennas at the BS, with 2 uplink users and 4 downlink 
users. The uplink and downlink path loss parameters are 
gf = -50dB,g^‘ = -A5dB, and gf = -56dB,gf = 
—(yldB^gf" = —65dB,gf‘ = —58dB. The inter-node inter¬ 
ference channel gain is ^ = —59dB, g{ 2 = —GOdB, 1 = 
—62dB,g2 2 = —55dB. Fig |2|a) plots the evolution of 
the utility function for proportional fairness utilities, i.e., 
Uk{r) = log r,k G U/C^j. Fig |2^b) plots the evolution 
of the utility function for minimum potential delay utility, i.e., 
= — l/r, fc £ U/Cd/- It can be observed that the 
algorithm converges at around 200 iterations. 


IX. Conclusion 

In this paper, we study a massive MIMO full-duplex system 
where the BS contains a large number of full-duplex-capable 
antennas serving multiple single-antenna half-duplex up and - 
downlink users. In such system, inter-node inference becomes 
the main design challenge, where transmission of uplink users 
creates interference to the reception of downlink users. Be¬ 
cause of the vast amount of inter-node interference channels, 
centralized control by aggregating all network information 
at the BS will incur significant overhead in the network. 
We develop an efficient and scalable one-hop information 
architecture that leverages the unique structure of the massive 
MIMO full-duplex system. Based on the one-hop informa¬ 
tion architecture, we propose a distributed power allocation 
algorithm to optimally manage the inter-node interference, and 
show how the information can be obtained via a overhearing- 
based scheme. We characterize the performance of the optimal 
power control algorithm in the asymptotic regimes where the 
number of users and the number of antennas at the BS all 
scale up. We show that power control in large number of 
BS antennas regime can not only bring utility gains, it also 
improves the energy efficiency of mobile users. 
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Appendix A 
Proof of Theorem 2 

We let \ Hence for all i € > 

We also let U'{r) = max,„{(C/™)'(r-), C/™ e U}. 

We prove this lemma by contradiction. Suppose there exist 
a subsequence lm,fn = 1,2 ,-and 0 < n < 1 such that 
]lUYm^ao\^f\/tCul = K- We let power 

allocation vector where = po-Pmox)* ^ with 

Po < p , and users 

i G and m being large, the achievable transmission rate 


when Pp = p^’-^^rn,new expressed as follows, 

' Mi^ ■ poP^Lgf ] 

No ) 

Mu ■ 




: 10 g( 




No 


-j + log( 


Mi„ 


ul ' 

^ i Jl 




(35) 


where gi{l) > po- From dTSl l and concavity of utility function, 
we have for i G 




>Ul{Rf{Pf-^-,MiJ)-\ogg,il,n). 


(36) 


For the downlink users, the achievable transmission rate 
satisfies 


pdl -pul ,'r< 


7 = iog( 


. pdljm,* dl 


Nui -nul,new 


I p^ 


No 


M ■ Pdl.l^-*gdl 


PijPoPmax + 


gtj + No 


J 

where 


(37) 




E 


iG/c: 


^ 91, 


gijPoPmax + Eiee‘m Pf’’"^'*9ij + -^o) 


=1 


gl.pr^* - gljpoP±x - Ni 


'Ul 


dijPO^rnax 

which satisfies. 


p 




* I 

gij 


■ No 


>1 + 


= 1 + 


(Eie-J.'- 9ij [pi - Po] P^ax) - No 

PijPPmax + Pf’^’^’* + Nq 

[p - Po] Pmax 9i3 - No 

PioPPmax + Eu 0 l pf'‘^'*9ij) + No 


Hence 


limsup/j (/m) >1 +limsup 

m—^oo m—¥oo 


[p-po]P:}!axE,e^^^ gij/Kil 

9l3pP^ax)/\F'ul\ 


[p-po]C^ 
- PE[4] 


where C^. = < a] with a being such that P{gij < 

a) = K. Next consider the difference in utility. From (lT7l) and 
concavity of utility, we have for downlink user j. 


- P,(Pf (Pf 






(38) 


We let pdi,i,n,new* (jgjjote the optimal downlink power 
allocation to maximize the sum-downlink-utility given uplink 
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power allocation p'u-i.im.new^ difference in overall 


utility is 

[Ht/f 

i,neWj^ ^ ^ 

1 

M 

’-■*+ Y, u'r* 

that lim sup 

lim inf | 

*e<T 



3&ic2t 

m—vco L 


1 ,new ^ ^ jjl-m ,new 

[Ef'b' 

*+ E uY* 

+ 5 



ie<T 

nic2t 


= [ 11 ^'.'- 

i,new jjlm,* _|_ 1 ^ ^ ^ 

jj-lm ,new 

jjlm,* 

> lim inf | 

m—¥co L 


According to the above lemma, we let wq and e be such 

I 

^ = e. Therefore 




‘’,Mi^))\ogpo 

-log/,(/„) 




Po 


ie/c: 






ul ( rtul^lm^fT'ew 


^im)) ■ ^ogg^{lr, 


1 , Ccuogf 

+ )-log/j(U 


--°'-^^^C^E[gi^]po' 


ie$ 


where the last inequality holds from (l36l l (l38T l. We then have 

liminf[ ^ Ui{Ki\pf'^-'^^'^,MiJ)-\ogg{lm) 

m—¥oo L A—' 

+ C/'«((Pf’'"‘’*,P“'’'™’”""'))-l0g/,Gm) 


=kIy liminf M,^)) logpo 

m—¥oo L 


dl 


°|/C'j| ^^C,E[gL]po 


■ log/jGm) 


> 0 , 


3^K.2 


> liminf 1 U'{Rf ,MiJ) logpo 

m—¥oo L ^ 




where the last equality holds because the first term in the 
limit decays to 0 with Im of ( |39] |, U'{r), U'{r) monotonically 
decrease with r, the path loss and because from the 

theorem statement liminfm_>oo \jiL > 0 . 

\^ul I 

Therefore choosing ^pdi,i,new*^^ asymptotically 

achieves utility higher than the optimal power allocation 
^■pui,i,* which contradicts to (p“ffbeing 
the optimal power allocation. We hence proved the Theorem. 


where recall that U'{r) = max^ U[{r). Note that 
lim sup Rf pul,Im,new ) 

m—¥cio 

l\^ pdl,lm,* dl 
IVI • i- • (JA 


< 


< 


m^oo 9ijPQP^ax + 

1^ . pdl,lm,* dl 
IVI ' -L j Uj 


lim sup log ( 

m—^oc ^ 

lim sup log ( 

m—>oc ^ 

lim sup log ( 

m —>( x > ^ 

lim sup log ( 

m—voo 


W n I nr. pul J-'p pul,lm,* I 

rnax ' i Uij 

111 


■No 


M-Pf’^^’*gf/\]Cl- 


dijPoPmax + No)/\K.. 


M-Pf'^^’*gfl\K};^i 

CuE[gUpQP^r,x 
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where recall that C was defined in the statement of the 
theorem, and was previously defined in this proof. We 
proceed with the following lemma. 

Lemma 1. Let = {j G /C^; : Pf''-’* < w^}/or w > 0. 
There exists w > 0 such that limsup;_^oo > 0. 

Proof: Suppose the statement is not true. Then for any 
w > 1 we have limi_>oo tfH- = 0. Therefore 

^ E pf’"'* ^ ^Ptot > Ptif 


l—foo 




Hence the total downlink transmission power exceeds the 
power constraint, thus establishing a contradiction. ■ 























